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ABSTRACT 

We present the first Spitzer Infrared Spectrograph observations of the 15.2 
/im bending mode of CO2 ice towards field stars behind a quiescent dark cloud. 
C0 2 ice is detected towards 2 field stars (Elias 16, Elias 3) and a single protostar 
(HL Tau) with an abundance of ~ 15 — 20% relative to water ice. C0 2 ice is 
not detected towards the source with lowest extinction in our sample, Tamura 17 
(Ay = 3.9 m ). A comparison of the Elias 16 spectrum with laboratory data 
demonstrates that the majority of CO2 ice is embedded in a polar H 2 0-rich ice 
component, with ~15% of C0 2 residing in an apolar H 2 0-poor mantle. This 
is the first detection of apolar CO2 towards a field star. We find that the CO2 
extinction threshold is Ay = 4 m ± l m , comparable to the threshold for water 
ice, but significantly less than the threshold for CO ice, the likely precursor of 
CO2. Our results confirm CO2 ice forms in tandem with H2O ice along quiescent 
lines of sight. This argues for CO2 ice formation via a mechanism similar to 
that responsible for H 2 ice formation, viz. simple catalytic reactions on grain 
surfaces. 

Subject headings: ISM: Lines and Bands, ISM: Molecules, astrobiology, astro- 
chemistry 
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1. Introduction 

Observations by the Infrared Space Observatory (ISO) demonstrated that CO2 ice is 
a ubiquitous component of the interstellar medium (ISM), with typical abundance of ~15 
- 25% relative to H 2 0, the dominant ice component (Gerakines et al. 1999; Nummelin et 
al. 2001). However, the origin of this common grain mantle constituent remains uncertain. 
Gas-phase production, with subsequent freeze-out, is believed incapable of reproducing the 
observed abundance. In the laboratory C0 2 forms quite readily via ultraviolet (UV) photol- 
ysis or ion bombardment of astrophysical ice mixtures; along with the thermal processing of 
ices, these mechanisms have been suggested as possible formation routes (d'Hendecourt et 
al. 1986; Sandford et al. 1988; Palumbo et al. 1998). Even if UV photons and cosmic rays 
are not present with sufficient fluxes, an alternative mechanism is grain surface chemistry 
(Roser et al. 2001; Frost, Sharkey, & Smith 1991, and references therein). Observations 
of grain mantles indicate that there is competition between hydrogenation and oxidation 
of atoms on grain surfaces (van Dishoeck 2004, and references therein). Nowhere is this 
more clear than for carbon-bearing molecules. C0 2 ice is less volatile than CO (Sandford & 
Allamandola 1990) and its formation locks carbon in saturated form on grains, stopping any 
subsequent hydrogenation towards more complex molecules on grains or in the gas. Thus, 
an understanding of the formation mechanisms of C0 2 ice is needed in order to understand 
the limits of grain surface chemistry in producing more complex organics. 

One way to constrain the origin of C0 2 ice in the ISM is to use bright field stars located 
behind molecular cloud material as candles that probe material remote from embedded 
sources, where, at suitably high extinction, ices are unlikely to be exposed to significant UV 
radiation or heating. ISO detected the C0 2 4.27 /zm z/ 3 stretching mode towards only two K 
giants, Elias 13 and 16 (Whittet et al. 1998; Nummelin et al. 2001). However, the detection 
of C0 2 ice towards any field star demonstrated that radiative processing of ices is unlikely 
to be responsible for the C0 2 ice production in quiescent material (Whittet et al. 1998). 

In this paper we report, and discuss the implications of, the observation of the z/ 2 bending 
mode of C0 2 ice at 15.2 /im towards 3 field stars and 1 protostar in the Taurus molecular 
cloud using NASA's Spitzer Space Telescope. 

2. Observations and Results 

Observations of each source (see Table 1) were obtained in 2004 using the Spitzer Infrared 
Spectrograph, hereafter IRS (Houck et al. 2004). Each object was observed in staring mode in 
the short-wavelength high-resolution mode spectrograph (Short-Hi; SH) which has coverage 
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from 10 to 19 fim with A/AA ~ 600. Each exposure was taken using a per cycle integration 
time of 120 seconds with 4 cycles. All data were reduced using the IRS Team's SMART 
program (Higdon et al. 2004), starting with data products from pipeline version 13. Our 
data reduction process mirrored that described by Watson et al. (2004) (see §2 in that paper). 

Figure 1 shows the IRS spectra with a clear C0 2 absorption feature at ~15 /im for 
Elias 16, Elias 3, and HL Tau. For Tamura 17, the star tracing the lowest extinction in our 
sample (Table 1), the stellar continuum is detected, but no absorption feature is observed. 
For CO2 the 15 /im feature optical depth is derived by fitting a multi-order polynomial to 
determine the continuum. In general polynomials with orders 2-4 (higher order for Elias 16) 
were required to account for the mismatch in the continuum on the longward side of the CO2 
line caused by wings of the broad 18.5 /im silicate feature. Column densities are estimated 
using N(C0 2 ) = / t(u ) dv / A 15 2f mi, where v is the frequency in wavenumber units and 
a value of 10~ 17 cm/molecule is assumed for A 15 _ 2 nm (Gerakines et al. 1995). The column 
density derived towards Elias 16 is identical, within errors, to that derived earlier by Whittet 
et al. (1998) using the ISO detection of the 4.27 /im. line. 



3. Elias 16 Line Profile 

The Elias 16 spectrum possesses a sufficient signal to noise ratio to permit an analysis 
of the ice composition (Gerakines et al. 1999; Nummelin et al. 2001). The 4.27 /iin feature 
detected by ISO constrained the majority of CO2 to lie within a single polar ice component 
with H 2 0:C0 2 :CO (100:20:3) at 20 K (Whittet et al. 1998). This contrasts with CO ice, 
which resides primarily in an apolar component (Chiar et al. 1995). 

The C0 2 bending mode is more sensitive than the stretching mode to the ice composition 
(Ehrenfreund et al. 1996), and the spectrum, shown in Figure 2, exhibits clear asymmetry. 
We have used a similar set of laboratory interstellar ice analogs (Ehrenfreund et al. 1996, 
1997) with a x 2 minimization routine (Gerakines et al 1999) to constrain the ice composition. 
Because the gas temperature of the Taurus cloud is ~10 K we only present fits using ice 
analogs at 10 K and have excluded any spectra that include appreciable amounts of 2 . 
Molecular oxygen has yet to be detected in the ISM either by direct detection in the gas- 
phase or by indirect methods on grain surfaces (Vandenbussche et al. 1999; Goldsmith et al. 
2000; Pagani et al. 2003). 

Our best fit (Fig. 2) requires two components, a broad feature consistent with C0 2 ice 
within the polar water mantle (H 2 0:C0 2 100:14) and a narrow component of C0 2 embedded 
in an apolar matrix (CO:C0 2 100:26). The majority of C0 2 resides in the polar component, 
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while apolar C0 2 ice accounts for ~15% of the C0 2 ice column. 1 

4. Ice Extinction Threshold 

It is known from previous studies of ice features in Taurus that the absorption strength 
correlates with extinction. The correlation line intercepts the extinction axis at a positive 
value, i.e. there exists a threshold extinction below which the ice feature is not seen, presum- 
ably because the grains in the more diffuse outer layers of the cloud are not mantled (e.g. 
Chiar et al. 1995, Whittet et al. 2001 and references therein). Figure 3 compares plots of 
column density vs. extinction for CO2 and CO. In the case of CO2, we combine both Spitzer 
and ISO observations. The field-star data suggest a correlation yielding a threshold extinc- 
tion Ay = 4 m ± l m , i.e. not significantly different from the value of 3.2" 1 ± 0.1 m reported for 
water-ice (Whittet et al. 2001). In contrast, the threshold estimated for CO (6.8 m ± 1.6 m ) 
appears to be significantly larger. These results are consistent with a model in which most 
of the C0 2 is in the polar H 2 0-rich component, whereas most of the CO is in the apolar, 
H 2 0-poor component. A larger threshold is expected for the latter because of its greater 
volatility, requiring a greater degree of screening from the external radiation field. 2 

Figure 4 plots ice-phase column densities for C0 2 vs. H 2 and CO. In the case of C0 2 
vs. H 2 0, there is a general trend - linear least-squares fits to field stars and YSOs are similar 
and intercept close to the origin. In contrast, the C0 2 vs. CO plot shows a tendency to 
divide into two distinct trends (Gerakines et al. 1999). For a given CO ice column, field 
stars have a lower C0 2 ice column than massive YSOs. Grains in front of field stars are 
covered by polar and apolar mantles. In contrast, in the warm envelopes of massive YSOs, 
the dominant factor is likely to be sublimation of apolar CO-rich ices (although some C0 2 
might also be produced by energetic processes). The lower-mass YSOs show an intermediate 
distribution. 



1 This apolar component is effectively hidden within the ISO 4.27 /zm spectrum of Elias 16 as there is 
little difference in the fits obtained with a purely polar mixture (Whittet et al. 1998) and those including a 
weak apolar component such as that found here. 

2 Note that the pre-main-sequence star HL Tauri does not follow the field-star trend in either frame of 
Fig. 3. Much of the extinction toward this object evidently arises in a circumstellar disk (e.g. Close et al. 
1997). Temperatures in the disk likely range from <~ 100K, where CO is entirely in the gas phase (see Gibb 
ct al. 2004a) to much higher temperatures where all ice mantles are sublimed. That CO2 is detectable in 
solid form toward HL Tau is consistent with its residence in a polar matrix. 
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5. Implications for CO2 Ice Formation 

The Spitzer data on field stars revealed that (1) most of the CO2 ice is embedded within 
the water ice mantle and (2) the CO2 extinction threshold is closer to the threshold for water 
ice than that of its presumed precursor molecule, CO. These two results strongly suggest 
that C0 2 ice formation occurs in tandem with that of water ice. Water ice is believed to 
form via surface reactions during phases when gas is rich in atomic hydrogen and atomic 
oxygen. Observations of water vapor imply that atomic oxygen must be depleted in dense 
evolved well shielded molecular regions; the low abundance of gas-phase water inferred by 
SWAS (Snell et al. 2000) and ODIN (Olofsson et al. 2003) can only be accounted for in 
models where nearly all available gas-phase atomic oxygen becomes locked on grains (Bergin 
et al. 2000). 

Thus H2O and CO2 ice formation must occur during the early lower density formative 
stages of the cloud. Recent models of molecular cloud formation behind shock waves by 
Bergin et al. (2004) may therefore be useful in setting constraints on ice formation. They 
found that H 2 formation occurs at earlier times than gas-phase CO formation because H 2 
efficiently self-shields while CO formation requires dust shielding (Ay ~ 0.7 mag). At such 
low cloud depths ice mantle formation would be retarded by UV photodesorption. However, 
for Ay > 1.0™ the effects of photodesorption are greatly reduced. In this scenario CO 
gas-phase formation precedes both CO and water ice mantle formation. This qualitatively 
answers several questions. Because most available gas-phase carbon would be locked in gas- 
phase CO, it would preclude a high abundance of methane ice, in accord with observations 
(Gibb et al. 2004b). It also allows for H, O, and CO to be present on the grain surface to 
react via simple catalytic reactions to create H 2 and C0 2 . 3 

One key question remains: how to account for the presence of apolar CO2 ice? If CO2 
forms in tandem with H2O in oxygen rich gas then how is a separate component of CO2 
formed with little H 2 0? This implies the presence of O I in gas with little H I. There are 
at least 2 scenarios that could account for this mantle structure: line of sight structure in 
the abundance of atomic hydrogen or atomic oxygen (or perhaps both). To examine the 
question of line of sight structure in atomic hydrogen, in Elias 16 the abundance of CO2 in 
the apolar mantle is ~ 3 x 10~ 6 (all abundances relative to H 2 ): thus gas-phase H I would 
need to fall below this value to stop O I hydrogenation. Atomic hydrogen is expected to 



3 CH30H ice is not detected towards field stars (Chiar et al. 1996). Thus observations would suggest that 
grain surface formation of CH3OH is inefficient under low temperature quiescent conditions. Laboratory 
experiments investigating the hydrogenation sequence for CO, CO — > HCO — ► H 2 CO — > CH 3 OH (excluding 
some intermediate products), are discussed by Hiraoka et al. (2005, and references therein). 
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have near constant space density in molecular clouds, nni ~ 1 — 5 cm -3 (see Goldsmith & Li 
2005). Thus, the H I abundance should inversely follow density variations along the Elias 16 
line of sight. However, the density would need to be > 10 6 cm -3 for the abundance of H I to 
fall below that required for O I. This density is characteristic of a condensed molecular core, 
which is not detected towards this line of sight (Cernicharo & Guelin 1987), and is therefore 
implausibly high. Thus H I abundance structure is insufficient to account for apolar CO2. 

Line of sight structure in the O I abundance provides a more plausible solution to this 
issue. If atomic oxygen were absent in the densest regions with high extinction, by water 
ice formation and other additional solid-state reservoirs, then oxygen hydrogenation could 
have halted in these regions. Oxidation could continue in layers with lower extinction and 
density. For instance, assuming a density of 10 4 cm -3 for the Elias 16 line of sight (Bergin 
et al. 1995), the atomic oxygen abundance would need to be > 10" 4 to be higher than H I, 
which is conceivable given the available O I (~ 8 x 10~ 4 ; Jensen et al. 2005). Thus CO 
oxidation could continue in outer layers rich in atomic oxygen. There is some evidence for 
large O I columns towards molecular clouds that may trace low density layers (Caux et al. 
1999; Lis et al. 2001; Li et al. 2002, and references therein). This qualitative model can 
be tested by future higher signal-to-noise Spitzer observations of low extinction field stars 
(e.g. Elias 3). For instance, line of sight structure in the oxygen abundance would predict 
the existence of apolar C0 2 at moderate optical depths, even at those below the CO ice 
threshold. Alternately, under the assumption that lower Ay implies lower density, then 
decreasing amounts of apolar C0 2 would suggest a relation to the declining abundance of 
atomic hydrogen with increasing density. 
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Fig. I. — Spitzer IRS Short-HI spectra towards background stars Tamura 17, Elias 3, and 
Elias 16; HL Tau is a Class I pre-main sequence star. These spectra are shown without any 
baseline subtraction. 
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Fig. 2. — Optical depth of C0 2 ice absorption seen towards Elias 16 (histogram). Best fit 
line profile using laboratory analogs of interstellar ice mixtures at 10 K (solid line). The 
best fit requires two components: C0 2 embedded in a water rich (polar) mantle, shown 
as the short-dashed line (laboratory ice analog mixture H 2 + C0 2 100:14). The second 
component, shown as long-dashed line is C0 2 embedded in an apolar mantle (laboratory ice 
analog mixture CO + C0 2 100:26). 
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Fig. 3. — Plots of ice column density vs. visual extinction for (a) CO2 and (b) CO. Open and 
filled symbols in (a) represent stretching and bending vibrational modes of CO2, observed 
with the ISO SWS (Gerakines et al. 1999; Nummelin et al. 2001) and SST IRS (this paper), 
respectively. The points for one object (Elias 16) observed in both vibrational modes are 
joined by a vertical dashed line. Triangles indicate upper limits. The circled cross on the 
Ay axis indicates the locus of the threshold extinction for H 2 0-ice (Ay = 3.2; Whittet et al. 
2001). The diagonal line in each frame is the linear least-squares fit to field stars. 
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Fig. 4.— Plots of iV(C0 2 ) vs. (a) iV(H 2 0) and (b) iV(CO) for all sources with data available 
from SST (this paper) or ISO SWS (Gerakines et al. 1999; Nummelin et al. 2001). The 
symbols have the following meanings: Taurus field stars (+); Galactic Center field stars (x); 
low-mass YSOs (open diamonds); high-mass YSOs (solid diamonds). The solid diagonal 
line in each frame is the linear least-squares fit to field stars. The dashed line in (a) is 
the fit to YSOs only (both high and low mass). The dashed line in (b) is the relation 
jV(C0 2 ) = 4.2 x 10 17 + l.lTV(CO) proposed by Gerakines et al. (1999) for high-mass YSOs. 
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Table 1. Observation Parameters and Results 



Source 


Obs. Date 


JV(£)Ai> a 


N(C0 2 ) b 


N(C0 2 ) 
N(H 2 0) 


N(C0 2 ) 
N(CO) 


Ay 


Elias 16 f 


Mar. 3, 2004 


5.0±0.1 


5.0 


0.21 


0.8 


23.3 m 


Elias 3 f 


Mar. 3, 2004 


1.8±0.1 


1.8 


0.20 


1.0 


10.0 m 


Tamura 17 f 


Feb. 27, 2004 


<0.1(3ct) 


<0.12 


<0.10 




3.9 m 


HL Tau 


Oct. 10, 2004 


1.9±0.1 


1.9 


0.14 




24.0 m 



a For Elias 16, Elias 3, and HL Tau integrated optical depths arc calculated by a 
direct integration over the profile. For Tamura 17, we have estimated the opacity 
limit by fitting a series of Gaussians with fixed width and line center determined 
by the Elias 3 feature. The absorption depth is variable and the minimum optical 
depth that fits the noise is used to estimate the 3a integrated opacity. Units are 
cm -1 . 

b In units of 10 17 cm" 2 . 

f Denotes field star. 

Note. — iV(CO) data for field stars and the limiting value for HL Tau are from 
Chiar et al. (1995) and Tegler ct al. (1995), respectively. The Ay value for 
HL Tau is from Close et al. (1997), those for field stars are calculated from the 
J — K color excess assuming Ay = b.'&Ej-K (Whittet et al. 2001). 



